Complex magnetic monopoles and geometric phases around diabolic and exceptional 

points 
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We study the geometric phase (GP)in presence of diabolic (DP) and exceptional (EP) points. 
While the GP associated with the DP is the flux of the Dirac monopole, the GP related to the EP, 
being complex one, is described by the flux of complex magnetic monopole. For open systems, in 
week-coupling limit, the leading environment-induced contribution to the real part of complex GP 
is given by a quadrupole term, and to its imaginary part by a dipolelike field. We find that the GP 
has a finite gap at the DP and infinite one at the EP. 
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For a quantum-mechanical system in the commonest 
case of double degeneracy with two linearly independent 
eigenvectors, the energy surfaces form the sheets of a dou- 
ble cone, and its apex is called a "diabolic point" (DP) 
In the context of Berry phase the DP is associated 
with 'fictitious magnetic monopole' located at the DP 
P, In opposition to the DP, at which the eigenval- 
ues coincide while the eigenvectors still remain distinct, 
exceptional points (EP) occur when the eigenvalues and 
eigenvectors coalesce. Since for Hermitian operators the 
coalescence of eigenvalues results in different eigenvec- 
tors, the EP degeneracy does not exists for the Hermitian 
operators and is associated with open quantum systems 
and non-Hermitian physics. 

While notion of the geometric phase (GP) for pure 
states is well defined, the definition of a GP in open 
quantum systems is still an unsolved problem. The first 
important step towards the consistent description of ge- 
ometric phase for open systems was done by Garrison 
and Wright They removed the restriction to unitary 
evolution considering the quantum system governed by 
non-Hermitian Hamiltonian. Uhlmann was the first to 
introduce the concept of GP to the case of mixed states 
employing the density matrix approach 0, Q . However 
the physical interpretation of the Uhlmann GP for den- 
sity matrix is not clear. Recently, an alternative defini- 
tion of the GP for mixed states has been proposed by 
Sjoqvist et al in the experimental context of quantum 
interferometry 

In this Letter we consider the GP in presence of DPs 
and EPs. We show that for general non-Hermitian sys- 
tem the related GP is described by complex magnetic 
monopole. We analyze GPs near the degeneracy and find 
that the GP is nonanalytical and, having a finite gap at 
the DP, diverges at the EP. 

General results on behavior of the eigenvectors at dia- 
bolic and exceptional points. - It is known that in param- 
eter space, a set of exceptional points defines a smooth 
surface of codimension 2 for symmetric/nonsymmetric 



complex matrix, codimension 1 for a real nonsymmet- 
ric matrix, and EP's do not exist for real symmetric or 
Hermitian matrix Q. 

Let H{X) be a complex N x N matrix smoothly de- 
pendent on m real parameters Xi (i runs from 1 to m). 
For Xk{X) being the eigenvalues of H{X), we denote 
by \ipk{X)) and (V'fc(^)| the corresponding right and 
left eigenvectors: -ffjV'fc) = >^k\i^k), {ipk\H = Xk{-ipk\- 
Both systems of left and right eigenvectors form a bi- 
orthogonal basis @]: 



\-ipk){ipk 



y 



1 {i>k\i^k) = o, k^k'. 



(1) 



Using the decomposition of unity ([T]), one obtains 
\^^)=E^aM and {^=E^P^{i'^l where 



(2) 



{ipi\il^i) {iPMi 

We assume that EP occurs for some value of param- 
eters X — Xc- At the EP the eigenvalues, say n and 
n + 1, coalesce: A„(Xc) — Xn+i{Xc), and the corre- 
sponding eigenvectors coincide yielding a single eigenvec- 
tor \ipEp) = \ipniXc)) = li^n+iiXc)). Now, applying (P) 
for fc = n and fc = n -|- 1 we find that at the EP the nor- 
malization condition is violated: (V'epIV'ep) = 0. This 
leads to the serious consequences for the global behavior 
of the states in parameter space. 

Since at the EP both eigenvalues and eigenvectors 
merge forming a Jordan block, it is convenient to in- 
troduce the orthonormal basis of the related invariant 
2-dimensional subspace as follows: 

(n|n) = l, (n+l|n+l) = l, (7i|n + l)=0 (3) 

Assuming that all other eigenstates are non-degenerate, 
we find that the set of vectors {(xfeL \ Xk)} , such that for 
k = n,n + 1 one has \xk) = \k), {xk\ = {k\, and 



\Xk) = 



{■>Pk\i'k 



{Xk\ = 



(ipkl^Pk) 



for k ^ n,n + 1, 
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forms the bi-orthonormal basis. Using this basis we ex- 
pand an arbitrary vector as \ip) = ^ Ck(X)\xk(X)) , with 
the coefficients of expansion being Cfc = (xfelV')- 

From the orthogonaHty condition, one can see that if 
\4'{X)) I'lpEp) while X Xc, then all (fc ^ n, n+1) 
vanish at EP. Thus in the neighborhood of EP only the 
terms related to the invariant subspace make substantial 
contributions and the A'^-dimensional problem becomes 
effectively two-dimensional 0, 0|- Similar conclusion is 
valid for the DP's, excepting that at the DP the eigen- 
vectors don't coalesce. The detailed study of the asso- 
ciated two-dimensional problem will be presented in the 
following sections. 

Geometric phase for non-Hermitian systems. - The GP 
for non-Hermitian systems were studied by various au- 
thors (see, for instance, [1 [IS III, [H, [H Q, [H and ref- 
erences therein). Following [3|, let us consider the time 
dependent Schrodinger equation and its adjoint equation: 



i^^\^(t)) = H{x{tm{t)), 

-i^^{^{t)\^{^{t)\H{X{t)). 



(4) 
(5) 



Let {tpn{X)\ and \^n{X)) be left (right) eigenstates 
corresponding to the eigenvalue i?„, then in adiabatic 
approximation the complex geometric phase is given by 
the integral 0, Q [H 



7« 



(^„(x)M^„(x))^ r 
Tc {MX)\MX)) Jc ^ ^ " 



Validity of the adiabatic approximation is defined by the 
following condition [SJ: 



E 



i-iPmldH/dtlMX)) 



{Em — En 



< 1 



(6) 



This restriction is violated nearby the degeneracies re- 
lated to any of DP or EP, where the eigenvalues coalesce. 

Since the adiabatic approach cannot be applied in the 
neighborhood of degeneracy, further we will consider non 
adiabatic generalization of Berry's phase introduced by 
Aharonov and Anandan [l^ and extended by Garrison 
and Wright Q to the non-Hermitian systems as follows. 
Let an adjoint pair {|\['(t)), (4'(t)|} being a solution of 
Eqs. (Ill, ^ satisfies the following condition: 

l^iT)) = exp{i^)\^iO)), ($(r)| = exp(-j^)(*(0)|, 

where (/? is a complex phase, and {|x(^))i (x(OI} ^ mod- 
ified adjoint pair such that 

\x{t))=eM~mmit)), (x(i)l -exp(*/(i))($(t)|, 

where f{t) is any function satisfying /(T) — /(O) = (p{0)- 
The total phase (p calculated for the time interval (0, T) 



may be written as ip ^ ^ + S, where the geometric phase 
7 is given by 

^-'[^X{t)\^^x{t))dt (7) 
and for the "dynamical phase" one has 

S ^ - r {m\H\xit))dt. (8) 
Jo 

This yields the connection one-form and the curvature 
two-form as follows: 



A = i{x\dx), F = dA. 



(9) 



Two-level system. - Let us consider a two- level system 
described by generic non-Hermitian Hamiltonian: 



H 



Xa + Z X-iY 
X + iY \a-Z 



(10) 



Further we consider a complex matrix El as depending 
on three complex parameters X, Y and Z. While the 
DP is just a point in 3-dimensional complex space C'^, a 
set of EP being a hypersurface of complex codimension 
1 in C^, defines a smooth surface of codimension 2 in 
6-dimensional real space 0, 01 ■ 

The solution of the eigenvalue problem 

H\u)^\\u), (u|i7 = A(u|, (11) 

where \u) and ({t| are the right and left eigenvectors, re- 
spectively, is given by A± = Aq ± ^ X'^ -f + Z'^ . The 
left and right eigenvectors are found to be 

-p-^V sin . 



cos ■ 



sm : 



(M_ = — e ^ sm — , cos — 
2 2 



,^^.+ 1= (e'^ cos-, sin-) (12) 



where 



cos - 
2 



R + Z 



2R ' ^^"2 



R- Z 
2R ' 



(13) 



^v ^ X + 'iY ^ + Y^ + Z2. (14) 

and 9, ip are the complex angles of the complex spherical 
coordinates {R,9,(p). For i? 7^ the following relation- 
ships hold: (u±|Mip) — 0, {u±\u±) = 1. 

As can be easily seen, the coupling of eigenvalues oc- 
curs when X'^ + Y'^ + Z'^ — 0. There are two cases: 
6 = 0, ip = and 9 = ioo, p = 0. The first case yields 
two linearly independent eigenvectors and the point of 
coupling is known as the diabolic point (DP). Since for 
9 ~f ioo the eigenvectors merge at the coupling point, the 
second case corresponds to the EP. At the DP we obtain 



l"+> 



1 


(1,0), 





1 

(0,1), 



(15) 
(16) 
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and at the EP 



(7i+| = (7i_|(x(e'^z).(17) 



This impUes the violation of the normahzation at the EP, 
and one has {u±\u±) = 0. 

Applying ([9]) to \u±) and {u±\, we obtain 

A± = =Fg(l±cos6')d(p, ^ q sine de /\ dip (18) 

where q = 1/2. This describes the complex "magnetic 
monopole" with a charge q and the field B given by 



B 



q 



R 

i?3 



(19) 



where R = {X, Y, Z) and R ^ \/X^ + Y'^ + Z^. As can 
be easily seen the field of the monopole can be written 
as = -d<l}/dX\ where the potential = q/R. 

Computation of geometric phase yields ^ — §^A, inte- 
gration being performed over the contour C on the com- 
plex sphere S^. Applying the Stokes theorem we obtain 



7 



/ F = qn{C) 



where E is a closed surface with the boundary C — 9S, 
and VL{C) is the complex solid angle subtended by the 
contour C. 

For DP formula (fT9|l reproduces the classical result by 
Berry on two- fold degeneracy in parameter space 0] . For 
EP the field of the corresponding "monopole" represents 
a complicated topological charge rather than a pointlike 
magnetic charge. 

Hyperbolic monopole.- Let us consider the following 
non-Hermitian Hamiltonian: 



H 



Xq + iz X — iy 
X + iy Xq — iz 



x^y.z^ 



The computation yields 
R 



B 



R = {x,y,iz) 



(21) 



(22) 



y2 _ z^, 



and B is the field of 



where R — ya 
"monopole" be called the hyperbolic monopole (FiglT|). 

The EP defined as the solution of the equation R — 0, 
is represented by the double cone with the apex at the 
origin of coordinates, and the DP is just located at the 
origin of coordinates. For R > we obtain the imaginary 
hyperbolic monopole {diB = 0): 

= =F'?(1 Tsinh0)dv3, F ^ qcoshO d0 A dip, (23) 

and for R < one has the real hyperbolic monopole 
(SB = 0)(see FigH]): 

A± = iFg(l ip cosh e)dip, F ^q sinh 6 d0 A dip. (24) 



FIG. 1: Imaginary hyperbolic monopole: the surfaces 3# = 
const are presented (left). Real hyperbolic monopole; the 
surfaces = const are plotted (right). 



Complex Dirac monopole. - For the complex Hamiltonian 



H 



Xq + z — ie X — iy 



■ le 



iy Xo- z 

the computation yields the field B as follows: 
qR 



x,y,zeR (25) 



B 



— , {x,y,z-ie) 



(26) 



(20) where R ^ yG 



y 



z^ — — 2iez. For e = the 



obtained complex Dirac monopole becomes well-known 
pointlike Dirac monopole. 

The EP being determined as the solution of equation 
x'^ + y"^ + z"^ — e"^ + 2iez = is the circle of the radius e a 
the plane z — 0. For e — I the real and imaginary parts 
of ^ = const are presented in Fig[2l 





FIG. 2: Complex Dirac monopole. The surfaces = const 
(left) and Q$ — const (right) are plotted. The EP is ap- 
peared as the circle of the radius r = 1 at the plane z — 



Let us write R as follows: R = \/r^ — 2ier cos x — £^ , 
where r — x'^ + + z^. We may expand 1/i? as 



\/r'^ ~ 2ier cos x ~ 



ESt^'(-«x)' (27) 



P/(cosx) being the Legendre polynomials. For r ^ e, 
this yields the following expansion of <?: 



^ ^ g ^ ^ pcosx _ Q(3cos^x - 1) 
r r^ 2r^ 



(28) 
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where q is the total charge of the monopole, p = qe is the 
dipole moment, and Q = qe^ the quadrupole moment. 

In the spherical coordinates (r, x, the GP for the 
state A_ = Aq — i? is given by 



1 ^ q J> I 1 



r cos X ~ 



\J — Her cos x — 



(29) 



Using multipole expansion 
lowing expression: 



7 = 7a/ + * 



psm X 



dip ■ 



, we obtain for 7 the fol- 
3Qsin^ Xcosx 



2r2 



where 7m = 'z/c^-'^ ~ cosx)'^'/' is the contribution of 
the Dirac pointlike monopole, second term describes the 
dipole contribution to the imaginary part of the geomet- 
ric phase and third term the quadrupole contribution to 
its real part. 

A complex Dirac monopole appears in wide class of 
open systems, where the Hamiltonian 



In summary, we show that, while the DP is associ- 
ated with Dirac magnetic monopole, the EP is related to 
the complex magnetic monopole. We found that for real 
part of GP the first correction to the flux of the Dirac 
monopole field is given by the quadrupole term, and the 
expansion for its imaginary part starts with the dipolelike 
field. Similar conclusion has been obtained for a two-level 
spin-half system in a slowly varying magnetic field and 
weakly coupled to a dissipative environment Note, 
our results are generic and should be applied to any dis- 
sipative system with an accidental energy levels crossing. 

Concerning experiment. - The singular behavior of the 
GP near the EP could be observed in the experiments to 
measure mixed-states GP with neutrons o with pola rized 
light, using a Mach-Zhender interferometer 2^, 2l[. 

The authors thank A. B. Klimov, J. L. Romero and 
S. G. Ovchinnikov for helpful discussions and com- 
ments. This work is supported by research grants SEP- 
PROMEP 103.5/04/1911 and CONACyT U45704-F. 



= B(i) • cr - -r^r 

includes spontaneous decay F = -^/ecr- as a source 
of decoherence. For instance, it emerges in a two- 
level atom driven by periodic electromagnetic field 
E{t) = 3fi(£(t) exp(iz/t)), with £{fy being slowly var- 
ied, as follows. In the rotatin g w ave approximation the 
Schrodinger equation reads 



lin g w a 



U2 



Ul 
U2 



(30) 



where 7^, are decay rates for upper and lower levels 
respectively, A ^ uj — i', lu = {Ea — -Eb), V — ^i* ■ £, and 
/X is the electric dipole moment. Removing the explicit 
time dependence of the Hamiltonian with the non unitary 
transformation [3] , one obtains H in the form of . To 
compare our results with that found in 0] we set x = 
W(t), y = 9^(0, z = A/2 and e = 5 = (7^ - 7b)/2. 
Let us consider the closed curve C parameterized by (f 
with the complex angle 6* = 6*0. Then the GP of Eq. (|29l) 
becoming 



7 



TT 1 



A-i5 



(31) 



V|2yoP + (A-z5)^ 

is the GP obtained by Garrison and Wright sf. 

Of the particular interest is the behavior of GP near 
the EP. For the resonance frequencies (A — 0) we obtain 



TT, if r > (5 

3*7 = <; , 5 



1± 



^) , if r < (5, 



(32) 



where the upper/lower sign corresponds to A ±0 and 
we set r = |2Vb|- This yields the singularity of GP at the 
EP {r ^ 5), and the GP has a finite gap at the DP [li|. 
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